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Shared Bottleneck Detection for Coupled Congestion Control for RTP Media 
Abstract 


This document describes a mechanism to detect whether end-to-end data 
flows share a common bottleneck. This mechanism relies on summary 
statistics that are calculated based on continuous measurements and 
used as input to a grouping algorithm that runs wherever the 
knowledge is needed. 


Status of This Memo 
This document is not an Internet Standards Track specification; it is 


published for examination, experimental implementation, and 
evaluation. 


This document defines an Experimental Protocol for the Internet 
community. This document is a product of the Internet Engineering 


Task Force (IETF). It represents the consensus of the IETF 
community. It has received public review and has been approved for 
publication by the Internet Engineering Steering Group (IESG). Not 


all documents approved by the IESG are candidates for any level of 
Internet Standard; see Section 2 of RFC 7841. 


Information about the current status of this document, any errata, 


and how to provide feedback on it may be obtained at 
https://www.rfc-editor.org/info/rfc8382. 
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1. Introduction 


In the Internet, it is not normally known whether flows (e.g., TCP 
connections or UDP data streams) traverse the same bottlenecks. Even 
flows that have the same sender and receiver may take different paths 
and may or may not share a bottleneck. Flows that share a bottleneck 
link usually compete with one another for their share of the 
capacity. This competition has the potential to increase packet loss 
and delays. This is especially relevant for interactive applications 
that communicate simultaneously with multiple peers (such as 
multi-party video). For RTP media applications such as RTCWEB, 
[RTP-COUPLED-CC] describes a scheme that combines the congestion 
controllers of flows in order to honor their priorities and avoid 
unnecessary packet loss as well as delay. This mechanism relies on 
some form of Shared Bottleneck Detection (SBD); here, a measurement- 
based SBD approach is described. 


1.1. The Basic Mechanism 


The mechanism groups flows that have similar statistical 
characteristics together. Section 3.3.1 describes a simple method 
for achieving this; however, a major part of this document is 
concerned with collecting suitable statistics for this purpose. 


1.2. The Signals 


The current Internet is unable to explicitly inform endpoints as to 
which flows share bottlenecks, so endpoints need to infer this from 
whatever information is available to them. The mechanism described 
here currently utilizes packet loss and packet delay but is not 
restricted to these. As Explicit Congestion Notification (ECN) 
becomes more prevalent, it too will become a valuable base signal 
that can be correlated to detect shared bottlenecks. 


1.2.1. Packet Loss 


Packet loss is often a relatively infrequent indication that a flow 

traverses a bottleneck. Therefore, on its own it is of limited use 

for SBD; however, it is a valuable supplementary measure when it is 

more prevalent (refer to [RFC7680], Section 2.5 for measuring packet 
loss). 
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1.2.2. Packet Delay 


End-to-end delay measurements include noise from every device along 
the path, in addition to the delay perturbation at the bottleneck 
device. The noise is often significantly increased if the round-trip 
time is used. The cleanest signal is obtained by using One-Way Delay 
(OWD) (refer to [RFC7679], Section 3 for a definition of OWD). 


Measuring absolute OWD is difficult, since it requires both the 
sender and receiver clocks to be synchronized. However, since the 
statistics being collected are relative to the mean OWD, a relative 
OWD measurement is sufficient. Clock skew is not usually significant 
over the time intervals used by this SBD mechanism (see [RFC6817], 
Appendix A.2 for a discussion on clock skew and OWD measurements). 
However, in circumstances where it is significant, Section 5.2 
outlines a way of adjusting the calculations to cater to it. 


Each packet arriving at the bottleneck buffer may experience very 
different queue lengths and, therefore, different waiting times. A 
single OWD sample does not, therefore, characterize the path well. 
However, multiple OWD measurements do reflect the distribution of 
delays experienced at the bottleneck. 


1.2.3. Path Lag 


Flows that share a common bottleneck may traverse different paths, 
and these paths will often have different base delays. This makes it 
difficult to correlate changes in delay or loss. This technique uses 
the long-term shape of the delay distribution as a base for 
comparison to counter this. 
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Definitions 
The key words "MUST", "MUST NOT", "REQUIRED", "SHALL", "SHALL NOT", 
"SHOULD", "SHOULD NOT", "RECOMMENDED", "NOT RECOMMENDED", "MAY", and 
"OPTIONAL" in this document are to be interpreted as described in 
BCP 14 [RFC2119] [RFC8174] when, and only when, they appear in all 
capitals, as shown here. 
Acronyms used in this document: 

OWD - One-Way Delay 

MAD - Mean Absolute Deviation 

SBD - Shared Bottleneck Detection 


Conventions used in this document : 


T the base time interval over which measurements 
are made 


N the number of base time, T, intervals used in some 
calculations 
M the number of base time, T, intervals used in some 


calculations, where M <= N 
sum(...) summation of terms of the variable in parentheses 


sum_T(...) summation of all the measurements of the variable in 
parentheses taken over the interval T 


sum_NT(...) summation of all measurements taken over the 
interval N*T 


sum_MT(...) summation of all measurements taken over the 
interval M*T 


ET) the expectation or mean of the measurements of the 
variable in parentheses over T 


E_N(...) the expectation or mean of the last N values of the 
variable in parentheses 


E_M(...) the expectation or mean of the last M values of the 
variable in parentheses 
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a) the count of measurements of the variable in 
parentheses taken in the interval T 


zat) the count of measurements of the variable in 
parentheses taken in the interval M*T 


a boolean variable indicating that the particular 
flow was identified transiting a bottleneck in the 


previous interval T (i.e., "Previously Bottleneck") 


t a measure of skewness in an OWD distribution 


skew_base_T a variable used as an intermediate step in 


var_est 


var_bas 


freq_es 


pkt_los 


calculating skew_est 


a measure of variability in OWD measurements 

e_T a variable used as an intermediate step in 
calculating var_est 

t a measure of low-frequency oscillation in the OWD 
measurements 

s a measure of the proportion of packets lost 


p_l, p_f, p_mad, c_s, c_h, p_s, p_d, p_v 


et al. 


various thresholds used in the mechanism 
number of values related to N 
ers and Their Effects 


T should be long enough so that there are enough packets 
received during T for a useful estimate of the short-term 
mean OWD and variation statistics. Making T too large can 
limit the efficacy of freq_est. It will also increase the 
response time of the mechanism. Making T too small will 
make the metrics noisier. 


N should be large enough to provide a stable estimate of 
oscillations in OWD. Often, M=N is just fine, though 
having M<N may be beneficial in certain circumstances. M*T 
needs to be long enough to provide stable estimates of 
skewness and MAD. 
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F determines the number of intervals over which statistics 
are considered to be equally weighted. When F=M, recent 
and older measurements are considered equal. Making F<M 
can increase the responsiveness of the SBD mechanism. If F 
is too small, statistics will be too noisy. 


c_s is the threshold in skew_est used for determining 
whether a flow is transiting a bottleneck or not. Lower 
values of c_s require bottlenecks to be more congested to 
be considered for grouping by the mechanism. c_s should be 
set within the range of +0.2 to -0.1 -- low enough so that 
lightly loaded paths do not give a false indication. 


p_l is the threshold in pkt_loss used for determining 
whether a flow is transiting a bottleneck or not. When 
pkt_loss is high, it becomes a better indicator of 
congestion than skew_est. 


c_h adds hysteresis to the bottleneck determination. It 
should be large enough to avoid constant switching in the 
determination but low enough to ensure that grouping is not 
attempted when there is no bottleneck and the delay and 
loss signals cannot be relied upon. 


p_v determines the sensitivity of freq_est to noise. 
Making it smaller will yield higher but noisier values for 
freq_est. Making it too large will render it ineffective 
for determining groups. 


Flows are separated when the 

skew_est |var_est|freq_est|pkt_loss measure is greater than 
p s|p mad|p flp d. Adjusting these is a compromise between 
false grouping of flows that do not share a bottleneck and 
false splitting of flows that do. Making them larger can 
help if the measures are very noisy, but reducing the noise 
in the statistical measures by adjusting T and N|M may be a 
better solution. 


2.2. Recommended Parameter Values 


[Hayes-LCN14] uses T=350ms and N=50. The other parameters have been 


tightened 


Section 4: 


F=20, and 
These are 
practical 


Hayes, et al. 


to reflect minor enhancements to the algorithm outlined in 
c_s=0.1, p_f=p_d=0.1, p_s=0.15, p_mad=0.1, p_v=0.7. M=30, 
c_h=0.3 are additional parameters defined in that document. 

values that seem to work well over a wide range of 

Internet conditions. 
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Mechanism 


The mechanism described in this document is based on the observation 
that when flows traverse a common bottleneck, each flow’s 
distribution of packet delay measurements has similar shape 
characteristics. These shape characteristics are described using 
three key summary statistics -- 


1. variability estimate (var_est; see Section 3.2.3) 
2. skewness estimate (skew_est; see Section 3.2.2) 
3. oscillation estimate (freq_est; see Section 3.2.4) 


-- with packet loss (pkt_loss; see Section 3.2.5) used as a 
supplementary statistic. 


Summary statistics help to address both the noise and the path lag 
problems by describing the general shape over a relatively long 
period of time. Each summary statistic portrays a "view" of the 
bottleneck link characteristics, and when used together, they provide 
a robust discrimination for grouping flows. An RIP media device may 
be both a sender and a receiver. SBD can be performed at either a 
sender or a receiver, or both. 


In Figure 1, there are two possible locations for shared bottleneck 
detection: the sender side and the receiver side. 


+=---+ 
| H2 | 
+=---+ 
| 
| L2 
+----+ L1 L3 +----+ 
| a1 |------ |------ EN 
4+----+ 4+----+ 


A network with three hosts (H1, H2, H3) and three links (L1, L2, L3) 
Figure 1 


1. Sender side: Consider a situation where host H1 sends media 
streams to hosts H2 and H3, and L1 is a shared bottleneck. H2 
and H3 measure the OWD and packet loss and periodically send 
either this raw data or the calculated summary statistics to H1 
every T. H1, having this knowledge, can determine the shared 
bottleneck and accordingly control the send rates. 
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2. Receiver side: Consider that H2 is also sending media to H3, and 
L3 is a shared bottleneck. If H3 sends summary statistics to H1 
and H2, neither H1 nor H2 alone obtains enough knowledge to 
detect this shared bottleneck; H3 can, however, determine it by 
combining the summary statistics related to H1 and H2, 
respectively. 


3.1. SBD Feedback Requirements 


There are three possible scenarios, each with different feedback 
requirements: 


1. Both summary statistic calculations and SBD are performed at 
senders only. When sender-based congestion control is 
implemented, this method is RECOMMENDED. 


2. Summary statistics are calculated on the receivers, and SBD is 
performed at the senders. 


3. Summary statistic calculations are performed on receivers, and 
SBD is performed at both senders and receivers (beyond the scope 
of this document, but allows cooperative detection of 
bottlenecks). 


All three possibilities are discussed for completeness in this 
document; however, it is expected that feedback will take the form of 
scenario 1 and operate in conjunction with sender-based congestion 
control mechanisms. 


3.1.1. Feedback When All the Logic Is Placed at the Sender 
Having the sender calculate the summary statistics and determine the 


shared bottlenecks based on them has the advantage of placing most of 
the functionality in one place -- the sender. 


For every packet, the sender requires accurate relative OWD 
measurements of adequate precision, along with an indication of lost 
packets (or the proportion of packets lost over an interval). A 
method to provide such measurement data with the RIP Control Protocol 
(RTCP) is described in [RICP-CC-FEEDBACK]. 


Sums, var_base_T, and skew_base_T are calculated incrementally as 
relative OWD measurements are determined from the feedback messages. 
When the mechanism has received sufficient measurements to cover the 
base time interval T for all flows, the summary statistics (see 
Section 3.2) are calculated for that T interval and flows are grouped 
(see Section 3.3.1). The exact timing of these calculations will 
depend on the frequency of the feedback message. 


Hayes, et al. Experimental [Page 10] 


RFC 8382 SBD for CCC for RTP Media June 2018 


Sis 


el 
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.2. Feedback When the Statistics Are Calculated at the Receiver and 
SBD Is Performed at the Sender 


This scenario minimizes feedback but requires receivers to send 
selected summary statistics at an agreed-upon regular interval. We 
envisage the following exchange of information to initialize the 
system: 


o An initialization message from the sender to the receiver will 
contain the following information: 


* A list of which key metrics should be collected and relayed 
back to the sender out of a possibly extensible set (pkt_loss, 
var_est, skew_est, and freq_est). The grouping algorithm 
described in this document requires all four of these metrics, 
and receivers MUST be able to provide them, but future 
algorithms may be able to exploit other metrics (e.g., metrics 
based on explicit network signals). 


* The values of T, N, and M, and the necessary resolution and 
precision of the relayed statistics. 


o A response message from the receiver acknowledges this message 
with a list of key metrics it supports (subset of the sender's 
list) and is able to relay back to the sender. 


This initialization exchange may be repeated to finalize the set of 
metrics that will be used. All agreed-upon metrics need to be 
supported by all receivers. It is also recommended that an 
identifier for the SBD algorithm version be included in the 
initialization message from the sender, so that potential advances in 
SBD technology can be easily deployed. For reference, the mechanism 
outlined in this document has the identifier "SBD=01". 


After initialization, the agreed-upon summary statistics are fed back 
to the sender (nominally every T). 


3. Feedback When Bottlenecks Can Be Determined at Both Senders and 
Receivers 


This type of mechanism is currently beyond the scope of the SBD 
algorithm described in this document. It is mentioned here to ensure 
that sender/receiver cooperative shared bottleneck determination 
mechanisms that are more advanced remain possible in the future. 


It is envisaged that such a mechanism would be initialized ina 
manner similar to that described in Section 3.1.2. 
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After initialization, both summary statistics and shared bottleneck 
determinations should be exchanged, nominally every T. 


3.2. Key Metrics and Their Calculation 


Measurements are calculated over a base interval (T) and summarized 
over N or M such intervals. All summary statistics can be calculated 
incrementally. 


3.2.1. Mean Delay 


The mean delay is not a useful signal for comparisons between flows, 
since flows may traverse quite different paths and clocks will not 
necessarily be synchronized. However, it is a base measure for the 
three summary statistics. The mean delay, E_T(OWD), is the average 
OWD measured over T. 


To facilitate the other calculations, the last N E_T(OWD) values will 
need to be stored in a cyclic buffer along with the moving average of 
E_T (OWD): 


mean_delay = E_M(E_T(OWD)) = sum_M(E_T(OWD)) / M 


where M <= N. Setting M to be less than N allows the mechanism to be 
more responsive to changes, but potentially at the expense of a 
higher error rate (see Section 4.1 for a discussion on improving the 
responsiveness of the mechanism). 


3.2.2. Skewness Estimate 


Skewness is difficult to calculate efficiently and accurately. 
Ideally, it should be calculated over the entire period (M*T) from 
the mean OWD over that period. However, this would require storing 
every delay measurement over the period. Instead, an estimate is 
made over M*T based on a calculation every T using the previous T’s 
calculation of mean_delay. 


The base for the skewness calculation is estimated using a counter 
initialized every T. It increments for OWD samples below the mean 
and decrements for OWD above the mean. So, for each OWD sample: 


if (OWD < mean_delay) skew_base_T++ 


if (OWD > mean_delay) skew_base_T-- 
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mean_delay does not include the mean of the current T interval to 
enable it to be calculated iteratively. 
skew_est = sum_MT(skew_base_T) / num_MT (OWD) 

where skew_est is a number between -1 and 1. 
Note: Care must be taken when implementing the comparisons to ensure 
that rounding does not bias skew_est. It is important that the mean 
is calculated with a higher precision than the samples. 

3.2.3. Variability Estimate 
Mean Absolute Deviation (MAD) is a robust variability measure that 
copes well with different send rates. It can be implemented in an 
online manner as follows: 

var_base_T = sum_T(|OWD - E_T(OWD) |) 
where 
|x| is the absolute value of x 


E_T(OWD) is the mean OWD calculated in the previous T 


var_est = MAD _MT = sum_MT(var_base_T) / num_MT(OWD) 


3.2.4. Oscillation Estimate 
An estimate of the low-frequency oscillation of the delay signal is 
calculated by counting and normalizing the significant mean, 


E_T(OWD), crossings of mean_delay: 


freq_est = number_of_crossings / N 


where we define a significant mean crossing as a crossing that 
extends p_v * var_est from mean_delay. In our experiments, we 
have found that p_v = 0.7 is a good value. 
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freq_est is a number between 0 and 1. freq_est can be approximated 
incrementally as follows: 


o With each new calculation of E_T(OWD), a decision is made as to 
whether this value of E_T(OWD) significantly crosses the current 
long-term mean, mean_delay, with respect to the previous 
significant mean crossing. 


o A cyclic buffer, last_N crossings, records a 1 if there is a 
significant mean crossing; otherwise, it records a 0. 


o The counter, number_of_crossings, is incremented when there is a 
significant mean crossing and decremented when a non-zero value is 
removed from the last_N_crossings. 


This approximation of freq_est was not used in [Hayes-LCN14], which 
calculated freq_est every T using the current E_N(E_T(OWD)). Our 
tests show that this approximation of freq_est yields results that 
are almost identical to when the full calculation is performed 
every T. 


3.2.5. Packet Loss 


The proportion of packets lost over the period NT is used as a 
supplementary measure: 


pkt_loss = sum_NT(lost packets) / sum_NT(total packets) 


Note: When pkt_loss is low, it is very variable; however, when 
pkt_loss is high, it becomes a stable measure for making grouping 
decisions. 


3.3. Flow Grouping 
3.3.1. Flow-Grouping Algorithm 


The following grouping algorithm is RECOMMENDED for the use of SBD 
with coupled congestion control for RTP media [RTP-COUPLED-CC] and is 
sufficient and efficient for small to moderate numbers of flows. For 
very large numbers of flows (e.g., hundreds), a more complex 
clustering algorithm may be substituted. 


Since no single metric is precise enough to group flows (due to 
noise), the algorithm uses multiple metrics. Each metric offers a 
different "view" of the bottleneck link characteristics, and used 
together they enable a more precise grouping of flows than would 
otherwise be possible. 
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Flows determined to be transiting a bottleneck are successively 
divided into groups based on freq_est, var_est, skew_est, and 
pkt_loss. 


The first step is to determine which flows are transiting a 
bottleneck. This is important, since if a flow is not transiting a 
bottleneck its delay-based metrics will not describe the bottleneck 
but will instead describe the "noise" from the rest of the path. 
Skewness, with the proportion of packet loss as a supplementary 
measure, is used to do this: 


1. Grouping will be performed on flows that are inferred to be 
traversing a bottleneck by: 


skew_est < c_s 
|| ( skew_est < c_h & PB) || pkt_loss > p_l 


The parameter c_s controls how sensitive the mechanism is in 


detecting a bottleneck. c_s = 0.0 was used in [Hayes-LCN14]. A 
value of c_s = 0.1 is a little more sensitive, and c_s = -0.1 is 
a little less sensitive. c_h controls the hysteresis on flows 


that were grouped as transiting a bottleneck the previous time. 
If the test result is TRUE, PB=TRUE; otherwise, PB=FALSE. 


These flows (i.e., flows transiting a bottleneck) are then 
progressively divided into groups based on the freq_est, var_est, and 
skew_est summary statistics. The process proceeds according to the 
following steps: 


2. Group flows whose difference in sorted freq_est is less than a 
threshold: 


diff(freq_est) < p_f 
3. Subdivide the groups obtained in step 2 by grouping flows whose 


difference in sorted E_M(var_est) (highest to lowest) is less 
than a threshold: 


diff(var_est) < (p_mad * var_est) 


The threshold, (p_mad * var_est), is with respect to the highest 
value in the difference. 
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4. Subdivide the groups obtained in step 3 by grouping flows whose 
difference in sorted skew_est is less than a threshold: 
diff(skew_est) < p_s 
5. When packet loss is high enough to be reliable (pkt_loss > p_l), 
subdivide the groups obtained in step 4 by grouping flows whose 
difference is less than a threshold: 


diff(pkt_loss) < (p_d * pkt_loss) 


The threshold, (p_d * pkt_loss), is with respect to the highest 
value in the difference. 
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This procedure involves sorting estimates from highest to lowest. It 
is simple to implement and is efficient for small numbers of flows 
(up to 10-20). Figure 2 illustrates this algorithm. 
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Figure 2 
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3.3.2. Using the Flow Group Signal 


Grouping decisions can be made every T from the second T; however, 
they will not attain their full design accuracy until after the 
2*Nth T interval. We recommend that grouping decisions not be made 
until 2*M T intervals. 


Network conditions, and even the congestion controllers, can cause 
bottlenecks to fluctuate. A coupled congestion controller MAY decide 
only to couple groups that remain stable, say grouped together 90% of 
the time, depending on its objectives. Recommendations concerning 
this are beyond the scope of this document and will be specific to 
the coupled congestion controller’s objectives. 


4. Enhancements to the Basic SBD Algorithm 


The SBD algorithm as specified in Section 3 was found to work well 
for a broad variety of conditions. The following enhancements to the 
basic mechanisms have been found to significantly improve the 
algorithm’s performance under some circumstances and SHOULD be 
implemented. These "tweaks" are described separately to keep the 
main description succinct. 


4.1. Reducing Lag and Improving Responsiveness 


This section describes how to improve the responsiveness of the basic 
algorithm. 


Measurement-based shared bottleneck detection makes decisions in the 
present based on what has been measured in the past. This means that 
there is always a lag in responding to changing conditions. This 
mechanism is based on summary statistics taken over (N*T) seconds. 
This mechanism can be made more responsive to changing conditions by: 


1. Reducing N and/or M, but at the expense of having metrics that 
are less accurate, and/or 


2. Exploiting the fact that measurements that are more recent are 
more valuable than older measurements and weighting them 
accordingly. 


Although measurements that are more recent are more valuable, older 
measurements are still needed to gain an accurate estimate of the 
distribution descriptor we are measuring. Unfortunately, the simple 
exponentially weighted moving average weights drop off too quickly 
for our requirements and have an infinite tail. A simple linearly 
declining weighted moving average also does not provide enough weight 
to the measurements that are most recent. We propose a piecewise 
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linear distribution of weights, such that the first section (samples 
1:F) is flat as in a simple moving average, and the second section 
(samples F+1:M) is linearly declining weights to the end of the 
averaging window. We choose integer weights; this allows incremental 
calculation without introducing rounding errors. 


4.1.1. Improving the Response of the Skewness Estimate 


The weighted moving average for skew_est, based on skew_est as 
defined in Section 3.2.2, can be calculated as follows: 


skew_est = ((M-F+1) *sum(skew_base_T(1:F) ) 
+ sum([(M-F):1].*skew_base_T(F+1:M))) 
/ ((M-F+1)*sum(numsampT (1:F) ) 
+ sum([(M-F) :1] .*numsampT (F+1:M))) 
where numsampT is an array of the number of OWD samples in each T 
(i.e., num_T(OWD)), and numsampT(1) is the most recent; 
skew_base_T(1) is the most recent calculation of skew_base_T; 1:F 
refers to the integer values 1 through to F, and [(M-F):1] refers to 
an array of the integer values (M-F) declining through to 1; and ".*" 


is the array scalar dot product operator. 


To calculate this weighted skew_est incrementally: 


Notation: F_ = flat portion, D_ = declining portion, 
W_ = weighted component 
Initialize: sum_skewbase = 0, F_skewbase = 0, W_D_skewbase = 0 


skewbase_hist = buffer of length M, initialized to 0 
numsampT = buffer of length M, initialized to 0 


Steps per iteration: 


T: old_skewbase = skewbase_hist (M) 
2. old_numsampT = numsampT (M) 

Shs cycle (skewbase_hist) 

4. cycle (numsampT) 

DE numsampT (1) = num_T (OWD) 
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6. skewbase_hist (1) = skew_base_T 
7. F_skewbase = F_skewbase + skew_base_T - skewbase_hist (F+1) 
8. W_D_skewbase = W_D_skewbase + (M-F) *skewbase_hist (F+1) 


- sum_skewbase 


oS W_D_numsamp = W_D_numsamp + (M-F) *numsampT(F+1) - sum_numsamp 
+ F_numsamp 


10. F_numsamp = F_numsamp + numsampT (1) - numsampT (F+1) 

11. sum_skewbase = sum_skewbase + skewbase_hist(F+1) - old_skewbase 
12. sum_numsamp = sum_numsamp + numsampT (1) - old_numsampT 

13. skew_est = ((M-F+1)*F_skewbase + W_D_skewbase) / 


( (M-F+1) *F_numsamp+W_D_numsamp) 


where cycle(...) refers to the operation on a cyclic buffer where the 
start of the buffer is now the next element in the buffer. 


4.1.2. Improving the Response of the Variability Estimate 


Similarly, the weighted moving average for var_est can be calculated 
as follows: 


var_est = ((M-F+1)*sum(var_base_T(1:F)) 
+ sum([(M-F) :1].*var_base_T(F+1:M) )) 
/ ((M-F+1)*sum(numsampT (1:F) ) 
+ sum([(M-F) :1] .*numsampT (F+1:M))) 


where numsampT is an array of the number of OWD samples in each T 
(i.e., num_T(OWD)), and numsampT(1) is the most recent; 
skew_base_T(1) is the most recent calculation of skew_base_T; 1:F 
refers to the integer values 1 through to F, and [(M-F):1] refers to 
an array of the integer values (M-F) declining through to 1; and ".*" 
is the array scalar dot product operator. When removing oscillation 
noise (see Section 4.2), this calculation must be adjusted to allow 
for invalid var_base_T records. 


var_est can be calculated incrementally in the same way as skew_est 


as shown in Section 4.1.1. However, note that the buffer numsampT is 
used for both calculations, so the operations on it should not be 
repeated. 
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4.2. Removing Oscillation Noise 


When a path has no bottleneck, var_est will be very small and the 
recorded significant mean crossings will be the result of path noise. 
Thus, up to N-1 meaningless mean crossings can be a source of error 
at the point where a link becomes a bottleneck and flows traversing 
it begin to be grouped. 


To remove this source of noise from freq_est: 


1. Set the current var_base_T = NaN (a value representing an invalid 
record, i.e., Not a Number) for flows that are deemed to not be 
transiting a bottleneck by the first grouping test that is based 
on skew_est (see Section 3.3.1). 


2. Then, var_est = sum_MT(var_base_T != NaN) / num_MT(OWD). 


3. For freq_est, only record a significant mean crossing if a given 
flow is deemed to be transiting a bottleneck. 


These three changes can help to remove the non-bottleneck noise from 
freq_est. 


5. Measuring OWD 


This section discusses the OWD measurements required for this 
algorithm to detect shared bottlenecks. 


The SBD mechanism described in this document relies on differences 
between OWD measurements to avoid the practical problems with 
measuring absolute OWD (see [Hayes-LCN14], Section III.C). Since all 
summary statistics are relative to the mean OWD and sender/receiver 
clock offsets should be approximately constant over the measurement 
periods, the offset is subtracted out in the calculation. 


5.1. Timestamp Resolution 


The SBD mechanism requires timing information precise enough to be 
able to make comparisons. As a rule of thumb, the time resolution 
should be less than one hundredth of a typical path’s range of 
delays. In general, the coarser the time resolution, the more care 
that needs to be taken to ensure that rounding errors do not bias the 
skewness calculation. Frequent timing information in millisecond 
resolution as described by [RICP-CC-FEEDBACK] should be sufficient 
for the sender to calculate relative OWD. 
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5.2. Clock Skew 


Generally, sender and receiver clock skew will be too small to cause 
significant errors in the estimators. skew_est and freq_est are the 
most sensitive to this type of noise due to their use of a mean OWD 
calculated over a longer interval. In circumstances where clock skew 
is high, basing skew_est only on the previous T’s mean and ignoring 
freq_est provide a noisier but reliable signal. 


A more sophisticated method is to estimate the effect the clock skew 
is having on the summary statistics and then adjust statistics 
accordingly. There are a number of techniques in the literature, 
including [Zhang-Infocom02]. 


6. Expected Feedback from Experiments 


The algorithm described in this memo has so far been evaluated using 
simulations and small-scale experiments. Real network tests using 
RTP Media Congestion Avoidance Techniques (RMCAT) congestion control 
algorithms will help confirm the default parameter choice. For 
example, the time interval T may need to be made longer if the packet 
rate is very low. Implementers and testers are invited to document 
their findings in an Internet-Draft. 


7. IANA Considerations 
This document has no IANA actions. 
8. Security Considerations 


The security considerations of RFC 3550 [RFC3550], RFC 4585 
[RFC4585], and RFC 5124 [RFC5124] are expected to apply. 


Non-authenticated RTCP packets carrying OWD measurements, shared 
bottleneck indications, and/or summary statistics could allow 
attackers to alter the bottleneck-sharing characteristics for private 
gain or disruption of other parties’ communication. When using SBD 
for coupled congestion control as described in [RTP-COUPLED-CC], the 
security considerations of [RTP-COUPLED-CC] apply. 
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